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Abstract 

Vte apply the MHD energy principle to the stability of a magnetized 
atmosphere which is boundef. below by much denser fluid, as is the solar 
corona. He treat the two fluids as ideal; the approximation which is 
consistent with the energy principle, and use the dynamical conditions that 
must hold at a fluid-fluid interface to Bhow that if vertical displacements of 
the lower boundary are permitted, then the lower atmosphere must be perturbed 
as well. However, displacements which do not perturb the coronal boundary can 
be properly treated as isolated perturbations of the corona alone. 
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I. INTRODUCTION „ 

t Studies of the equilibrium and stability of magnetized plasma in a 
gravitational field tire important to many areas of astrophysics, including 
•olar and stellar physics (reviewed by Priest, 1982, and Rosner et al,, 
1984). The stability of structures in the so'.lar corona is relevant to 
understanding the onset of eruptive activity, as well as the necessary 
conditions for equilibrium. Because even the simplest models of coronal 

features are sufficiently inhomogeneous that solving the full mode problem is 
very difficult, many studies of coronal MHD stability have used the energy 
principle method of Bernstein et al. (1958; hereafter BFKK) to determine 
stability without calculating the modes themselves. 

Coronal magnetic fieldlineo are thought to be connected to the lower 
solar atmosphere (chromosphere and photosphere) , which is much denser than the 
corona, and ultimately to extend into the solar interior. Rather than 
considering the stability of the composite system consisting of hot gas and 
cooler underlying material, most studies of coronal MHD stability have imposed 
a boundary at the coronal base and have treated the lower atmosphere only 
through its influence on the boundary conditions. 

Several different assimptions about the boundary conditions on £|, the 
component of the fluid displacement t parallel to the magnetic field have 
been made in the literature. Schindler et al. (1983) chose t = 0, as if the 
photosphere were a rigid boundary. Einaudi and van Hoven (1981) imposed 
parity constraints on £| that allow €| f 0. Hood 0 9B4a,b) did not explicitly 
restrict 1 \ at all. 

In this paper, we discuss the influence of the photospheric boundary 
condition on stability by assuming that both the upper and lower atmosphere 

are ideal fluids. It is clear that material in "the solar acnosphere does not 
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always behave adiabatically . Radiative processes, thermal conduction, and some 
fora of heatinf all play roles in the structure. Maos flow is often present. 
The upper and lower fluid treatment in this paper is an idealization, but it is 
an idealization which is consistent with the MHD energy principle, and should be 
a good approximation as long as the MHD time-scales are rapid compared to the 
timescale on which mass exchange occurs between the two fluids. Equilibria with 
flows, and non-adiabatic perturbations, cannot be studied with the ideal MHD 
energy principle. 

We find that for coronal stability problems in which gravitational 
stratification is included, the effect of nonzero 5| is to force the lower 
atmosphere to be perturbed. This arises in a natural way from the conditions at 
the boundary between the two fluids. Therefore, in order to derive necessary 
and sufficient conditions for the stability of the corona as an isolated system, 
must vanish on the lower boundary. 

In Section II, we use the MHD energy ]?rinciple to demonstrate the 
existence of a surface integral and a perturbation of the lower atmosphere 
when £| is not zero at the boundary. In Section III, we discuss the effect of 
the boundary terms on various results in the literature. Section IV is a 
discussion together with conclusions. 

II. BOUNDARY CONDITIONS AND ENERGY PRINCIPLE ANALYSIS 

We first describe the equilibrium model, including the conditions which 
■ust be fulfilled at the coronal base. We then derive the corresponding 
conditions in the presence of small perturbations. Finally, we us® the MHD 
energy principle of BFKK to assess the effect of the boundary conditions on 
stability. Ideal MHD (adiabatic, inviscid, infinite electrical conductivity) 
is assumed to hold throughout. 
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a) Boundary conditions 

The equation of mechanical equilibrium in a stratified atmosphere is 


V ■ [b B - l(p + B 2 /2 ) J + pg m 0 


( 1 ) 


where P, 6, P, and g are the gas pressure, magnetic field, gas density, and 
gravitation ill acceleration, respectively, Fbr a surface of discontinuity with 

A 

normal direction n in the fluid, it follows from V • b = 0 and from equation 
(1) that 


<B > « 0 
n 


(2a) 


<B B > » 0 
n t 


{(b 2 “ B 2 J/2 - p) - 0 


(2b) 


(2c) 


0 


where the notation <S> refers to the jump in quantity S across the 
discontinuity and the subscripts n and t refer to normal and tangent 
directions to the surface, respectively. Note that if B fi » 0, B t may be 

discontinuous, but if B n j 1 0, B n , B t , and P are each continuous separately. 
These conditions are discussed, e.g. by Roberts (1967). 

Now consider a small displacement t of the fluid. The Euler ian 
perturbations of P, S, and p are (BFKK) 


fip«-YPV»f-|»Vp 


6b - V x (£ x b) £ Q 


(3a) 


(3b) 


o 


o 
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(3c) 


Faraday's law implies that, to rirst order in £, 

<n * (t x b)> « o , , (4) 

Furthermore, t n must be continuous. 

Conditions for the perturbed system analogous to equations (2) were 
derived by BFKK for a fluid-vacuun interface tangent to £, and by Goedbloed 
(1979) for a fluid-fluid interface, again with B * n • 0, while Roberts (1967) 
gives conditions for % • n ji 0. Equations (2) must be linearized and satisfied 
at the perturbed boundary, with reference to the perturbed normal. The 
relevant perturbations of the fluid variables here are the Lagrangian 
perturbations, which follow the boundary elements to their new positions. The 
Lagrangian perturbations Ap, &,%, and Ap can be obtained from their Eulerian 
counterparts (3a,b,c) by the usual relationship for any quantity 

As « As + t • -Vs . (5) 

As Roberts (1967) shows in detail, for the case % • n / 0, 

<Ab> * <Ap> - o (6) 

Equation (6) is the analog of equations (2) above; the Lagrangian 
perturbations of P and % as veil as F and 3 themselves are continuous across 
the interface, 
b) Energy principle 



where the integral extends over the volume of the fluid. In the general case, 
both the photospheric and coronal fluids contribute to dw. 

BFKK proved (see also Freidberg 1982) that the system is unstable if and 
only if 6w(t,t) is negative for a displacement vector t which satisfies 
appropriate boundary conditions. In the most restricted sense, t must satisfy 
continuity conditions such as equations (4) and (6) and continuity of £ n at an 
interface. However, BFKK proved an extended energy principle for the plasma- 
vacuum problem with & • n = 0. They showed that t need not satisfy the 
Iagrangian force-balance condition (equation (2.32) in their paper; equation 
(6) here) by proving that it is possible to correct t in a thin layer near the 
interface in a way which enables t to satisfy the pressure balance 
condition. In their construction [see also Roberts (1967)), t is augmented by 
a vector cn which goes to zero within a distance e from the boundary. Then, 
the normal gradient of n is of order , and the contribution of eft to 6vfp is 


o 
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of order e. the extended energy principle makes it possible to choose trial 

functions t for 6w which do not satisfy the force balance condi tion. Roberts 

(1967) diacuBses the extended energy principle in cases where $ 1 n i 0. 

Since the Lagrangian force balance condition involves both t and itB 

derivatives, and takes some care to satisfy, the extended energy principle is 

easier to work with and is used in most applications. 

» 

When the extended energy principle iB written in its usual form (BFKX; 
Bq. 3.16), surface integrals involving the pressure balance condition 
appear. These integrals are related to the change in plasma potential energy 
caused by the PdV work done at its surface, we now consider the role of these 
boundary terms. According to equations (8) and (9), 


26w(t,t) - -Jd 3 x [t*v(ypv*t + e*Vp) + 6 *(Vxq)xb + t* ($*b)*q - t*g V'pif] 


( 10 ) 


Integrating by parts, this can be written 
- 26 w + 2dw 

£ S 

2foy= / d 3 x[(v»t) (yp + t*vp) + C 2 ~ I'(Vxb)xq + t*g V«pt] 

26w g - - J d 2 x ((n*t) (YP V*| + £-Vp) + Q • (n * (£*!)) ) OD 

The surface integral 6w s is taken over the boundary between the upper and 
lower atmosphere, plus terms at infinity, which we assume vanish. If the 
horizontal extent of the structure is finite, we can Impose horizontal 
periodic boundary conditions. The volume integral 6Wp contains contributions 




ui 


from noth fluids. 

Proceeding similarly to BFKK (see also Roberts, 1967) we rewrite the 
surface integral in equation (11) using the boundary conditions satisfied by 
t. Note that even if we use the extended energy principle, so that we allow 

trial functions in the volume integral fiv? F which do not satisfy the boundary 

* 

conditions (6), we must evaluate 6w s assuming that these conditions are 
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satisfied. This has not always been done in the literature. 

According to equation (4), n * (t * t) is continuous at the interface. 

Using equations (5) and (6), Ap and $ + l * are continuous as well, as 
♦ * 

is t • r». lhese results enable us to write 

2fiw g - - £ d 2 x((n«£) [£*V(p + B 2 /2 )] - (n»B) (e*Vb)'£) 

Using the equation of mechanical equilibrium, this becomes 

26w - £ d 2 x((n*£) [(b*Vb)* 5 + p£*g] - (n*B) (£*Vb)*<;). 
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or 


2«W S « - f> d 2 x((n*1){|*g)p+ [n*(£*B) »V b] ■£ ) 

Evidently, the second term involves only tangential derivatives of S at the 
interface. But the tangential derivatives of £ are continuous; this term is 
therefore zero. If we take g ■= yg and let th J . boundary lie in the x-z plane, 
then takes the final form 


26 w s « / dxdz^gCp^ - P u ) 


i, 


9 

I 

where p L and P u denote the densities in the lower and upper atmosphere, 
respectively. 

Equation 0 2) is exactly what is expected when one considers the 
Rayleigh-Taylor instal ility between two media of different density <e.g. 
Chandrasekhar, 1961). Since P t >> P u in the problem considered here, the 
surface term is positive. Thus we have shown that for nonzero g, the presence 
of flow across the unperturbed fluid boundary (nonzero £ t ) tends to be 
stabilizing. The surface term given in equation (12) arises naturally from 
the dynamics of the problem, and must be included in any evaluation of <5w. 

We can write 

6w - 6v? pc + 6w pL + 6w g 

where 6Wp C and fiWp L are the contribution of the corona and the lower 
atmosphere to Cw F as given in equation (11), and 6w s is the surface term given 
in equation (12). 

It is clear that if the problem of the stability of isolated coronal 
structures has any meaning we must be able to make 6Wp L vanishingly small. 

This requires in general, that t be nonzero only in an infinitesimally 
thin layer below the boundary, can the argument used in deriving the extended 
energy principle be applied to this situation? That is, is there always a 
displacement t of the lower atmosphere which satisfies the interface 
conditions but which makes 6 w FJj arbitrarily small in magnitude? 

In general, there io not. Recall from the discussion of the extended 
•nergy principle following equation (9) that the correction vector to 5 is 
assumed to be of order e and localized to a layer of width e. its 
contribution to 6Wp is then of order c. But in, the present case, if t *n is 
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of order 1 at the interface, and t IB again localised to a layer of order c, 
6 w FJj will be order s’* 1 '. We cannot always make 6w pij negligibly small. 

Hie case t • n ■ 0 is an exception. In this case, t can be zero on the 
boundary and of order c in a layer of width c. For example, take y « 0 to be 
the surface) y + -•* with depth. Take the iJigrangian pressure and magnetic 
field perturbations Ap s (x,z), aS s (x,£) to be proscribed by the displacement of 
the upper atmosphere, ttien, for y < 0, let 


5 x (x,y,t) 

S y (x,y,r) 

E z (x,y,s) 


f K, - ^r* 5 " W) 

y 


cAp 

"W 


1 e y/c (l-. y/c ) 


y 


where the functions multiplying the exponentials are evaluated at y ■ 0. This 
choice of t gives 6w p >of order e £c.r y < 0. 

Obese arguments have the following implications for MHD stability 
analysis of the solar corona. If we restrict ourselves to displacements 
* n ■ 0, then oW g and fiWp^ can both be nade zero, and it is both 
necessary and sufficient for the stability of the "isolated" coronal modes 
that $W pc > 0. If we consider displacements with t • n j< 0, the positive 
definite term Aw s (cf. equation 12) must be added to fiw pc . In addition, 6w pL 
must also be minimized. ThiB requires an explicit model of the lower 

atmosphere, but none of the presently available coronal equilibrium models 
include the lower atmosphere. OhuB, the stability of modes with t • n f 0 is 
Indeterminate. In summary, it is sufficient but not necessary for the stability 
of the isolated (t * n ■ 0) displacements that is positive when minimized 
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with unreotricted £ • n. Modes with £ • n * 0 cannot properly be tested for 
stability without a model of the lover atmosphere. 


III. EFFECT OF SURFACE TERMS ON CORONAL STABILITY 
In thiB section, we discuss the relevance of the term $w s ^ w fc 

derived in Section II to a number of studies of coronal MHD stability in the 

* 

literature in which different lower boundary conditions were assumed. We 

first consider the conditions originally used by Einaudi and Van Hoven (1961) 

and then discuss the conditions used by Schindler et al. (1983), Finally, we 

treat the work of Hood (I984a,b), for which the necessary analysis is somewhat 

•¥ 

wore involved, in all the papers we will treat, the components of £ 

perpendicular to the magnetic field axe assumed to vanish at the base of the 
atmosphere. The physical motivation for this is that the fieldlinos are 
assumed to be fixed in dense, infinitely conducting photospheric gas. 

a) Line Tying with Flow at Lower Boundary 

i 

Einaudi and Van Hoven (1981) studied the stability of coronal loops 
idealised as cylinders of finite length with twisted magnetic fields. They 
did not include gravitational stratification, an approximation which applies 
when the thermal scale height much exceeds the size of the system. the 
conditions they imposed at the ends of the cylinder are 

£_L ■ 0 at z ■* ±L 


£,(-L) - £| (L) 



£ <-L) 
£ 


I? «.«■> 


(13) 




• (a 

1 2 


* 

where r. is along the axis of the cylinder and Cj and £j_ are the components' of 
t parallel and perpendicular to B. These conditions have also been used by 
Migliuolo et al, (H’B4), Unaudi and Van Haven, (1 983) and references 

t * 

•n jt o in their model, they can derive sufficient but not 
necessary stability conditions for the isolated nodes. 


b) Rigid Boundary Condition 

Schindler et al. (1983) used the condition t m 0 on the lower boundary, 
which corresponds to treating the photosphere as a rigid, perfectly conducting 
wall. It is clear in this case that 6W S and $W FC vanish. Thus, none of the 
stability results arising from their minimisation of 6w F are affected by 

addition of the corresponding fiW s orftf FC . In a sense, the t “ 0 boundary 

♦♦ 

condition is a limiting case of the two fluid analysis for + ", If + ", 
fiw s becomes large unless £ n + 0. One must also impose “ °» 

c) Line Tying with Unrestricted Z \ 

Hood (1983a) derived a form of 5w including gravitational 
•tratification. He used the lower boundary condition ? x £ ** 0, although only 
in the approximate sense described below, and did not restrict £|. Although 
his formulation (as does that of Schindler et al.) extends to systems in which 
the magnetic fields have three spatial components, the analysis here is 
restricted to systems in which the fieldlines lie in parallel, ver‘ *1 

planes, such that 


£ , Z 3 * * 

B - + x i7“ y-^ 


(14) 


All the equilibrium quantities are functions of x and y only, and the magnetic 
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fieldlines are assumed to form loops which are symmetric in x. Ihese systems 
resemble solar magnetic arcades. Defining A^ a* t • ^A and using (14), Hood's 
boundary condition becomes Aj ■ 0 on the boundary. Using only this boundary 
condition, he writes $W FC as 


3a 


2% c - /A (-sr-) 2 ^ (»» ’ »(— j - «]» /+ (»•»«,) 2 * 


2jT 8J' 

7 


[B VA V V^ .,]2 + p[(y-i)(V . I} 2 + e 2y/H (V • te‘ y/H ) 2 ] (15) 


Here, J is the current density, H is the thermal scale height, 
and 3/3s ■ 1/B I • V is the derivative along a fieldline. 

It derive stability criteria, Hood assumes that the perturbations are 
isothermal (Y M 1) and minimixes with respect to 3£ z /3z and V •(£e'" y / H ). Ohis 
ref in the conditions 


V * t 



(16) 


3£ 

z 

rr 


( VVVA 1 «A 1) 


B 


(17) 


Equation (16) is a generalization of the incompressibility condition that 
results from minimizing 6W in plasmas that are not gravitationally stratified 
(e.g. Freidberg, 1982). Equations (16) and (17), taken together, imply that the 
total (gas plus magnetic) pressure perturbation is zero, and thus eliminates the 
stabilizing restoring force due to a pressure perturbation. 

Hood then takes the limit of infinitely large wavenumber in the z 
direction, k t ♦ Since condition (17) requires that the product k z £ z be 
finite, the magnetic tension term (B»V6 t ) 2 in 6W, which corresponds to 


^ WFYM**** '? ?-* **» .^»?J ** ^ * 

7 “’ " ■ 

bending the fieldlineo out of their equilibrium plane, becomes negligible, (see 

9 

Gilman 1970, Asseo et ul,1980, Zweibel 1981). Thus, 6W is reduced to the form 

2«W - /d 3 > r (('g~) 2 +[VA«V(£j)- i£ |~] A^ 2 }, 0 8) 

B B 

which is considerably simpler than equation 0 5), since only the perturbation 
variable A^ and its derivative along a fieldline appear. When 6w in the form 
(10) is minimised subject to the normalization condition /d 3 x Aj 2 “ 1, the 
resulting Euler equation for Ai is an eigenvalue equation. The solution of 
the eigenvalue problem is the basis of Hood's stability analysis (Hood 1984b) 
of the Zweibel and Rundhausen (1982) aquilibriun solutions. 

The minimizing conditions (16) and (17), together with the assumption that 
Y ** 1 and {^ ■ 0 at the lower boundary, guarantee that 6Ws as it appears in 
equation (11) vanishes. That is, Hood's displacements have a nonzero Lagrangian 
pressure perturbation on the lower boundary. As we Bhowed in Section II, the 
appropriate form of 6Wg is really Eq. (12), because of the Lagrangian force 
balance condition. Since €| is unrestricted, 6Wg will generally not vanish, 
and dWpi, cannot vanish either. 

We should also note that when 3£ s /3z is given by equation (17), 5 Z will 
not vanish on the lower boundary (because Va^ does not). Therefore, the 
condl tion Ej ■ 0 is technical! v violated. However, since £ z is arnall (Bee the 
discussion following equation (17)) it appears consistent to neglect £ z on the 
boundary when dropping (B • V£ z ) 2 from 6w pc . 

Suppose that when 6w p as given by equation (10) is minimized over Aj , the 

•f 

resulting £ (which can be calculated from A 1 using conditions (16) and (17), 

» * 

as we do below) aatisfies £ • n ■ 0. In this case, 6 w p q itself will be a true 
minimum. art the other hand, if t * n f 0, the minimization is unacceptable 
because 6w £ and 6w pL must be included. 
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t * 

We now discuss the conditions under which the miniroizati&n of equation 

. A 

(18) will permit T • n ■> 0. Since 6w FC in equation (18) only contains 
derivatives of A 1 along a fieldline, we can consider perturbations which are 
localized to a single flux-tube. To solve for t • n In terms of any given A^ , 
we eliminate d£ z /^z between conditions (16) and (17) to give 


(B • V?) • B - (B • Vb) 



(19) 


Then, E x may be written in terms of Aj and C y and eliminated from (19). The 
result is a first order differential equation for E y as a function of A 1 , with 
solution 


E (a) « B (s)e y/H — «' y V (— - - £*-J. 

y ' ' B y'* •* B l B0 ds B ds B 

7 c y y 


( 20 ) 


We assume here that B^ does not vanish anywhere on the field line, so s is a 
single valued function of xj the fieldlines Hood studied have these 
properties. Let the fieldline end tit ±s 0 . Then, if vanishes at the 
andpoints, equation (20) implies that 


f*0 djs|_ -y’/H fit d A i A i dB x i m 
' B '■B 3s B BB ds 

"*« y y 


t21 ) 


Since, in the geometry assumed, B x and By, are even and odd functions of x (or 
b) , respectively, equation (21) will be satisfied only if A^ is an odd 
function of x. 

We now consider Hood's study of the Zweibel and Hundhausen (1982) 
equilibria (Hood 1984b). These equilibria form a one parameter family in 
which the parameter measures the volume electric current, or distortion of the 
fieldlines from a potential field at the base of the atmosphere. The only 
previous stability analysis of these equilibria was a local analysis (Zweibel 

I 

1981) which showed that some portion of all the ZH equilibria were locally 
unstable. However, this analysis did not consider the stabilizing effect of 
magnetic, tension. Hood found that instability along an entire fieldline only j 



exiBts if the parameter 2<*H which measures the current exceeds a certain 
threshold. His solutions for A-) are even , rather than odd, functions of x. 
According to the arguments above, the stability boundary for iaolated coronal 
modes should be at a larger value of the current parameter than that found by 
Hood. 1 The stability boundary for these modes can be found by solving the 
Euler equation for the Integral (16) with Aj ■ 0 at the end and apex of a 
fieldline. Within the framework of the present analysis, Hood has found a 
sufficient, but not necessary, stability condition. 


. IV. DISCUSSION AND CONCLUSIONS 

In this paper, we have considered the .lower boundary condition for 
coronal MHD stability problems. lhese systems are characterized by magnetic 
fieldlines which connect the corona to a much denser, underlying atmosphere. 
Iheir ideal MHD stability has been studied using the BFKK energy principle. 

The lower atmosphere has simply been modeled as a rigid, conducting wall in 
some previous treatments (e.g. Schindler et al. 1983) on which the fluid 

4 • 

displacement £ vanishes. Other studies have allowed a non-vanishing fluid 
displacement parallel to the magnetic field at the lower boundary, but required 
the perpendicular components 5^ to vanish. 

In order to understand the effect of different boundary conditions, we 
considered the boundary as a contact surface between two ideal fluids of 
different temperatures and densities. We reviewed the boundary conditions which 
apply to a fluid-fluid interface and pointed out that these boundary conditions 
lead to surface Integrals in the perturbed potential energy 6W which represent 
PdV work done at the interface between the fluids. These terms (Eq. 12) can be 
written in a form which involves the density contrast between the two fluids, 
and is the same term one derives in an analysis of the Rayleigh Taylor 

1 After this paper was accepted for publication, we became aware of a recent 
study by Hood (1984c) which is consistent with this. 
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instability at an interface between unmagnetized fluids. Since the lower 
atmosphere is auch denser than the upper atmosphere, the surface term is 
strongly stabilizing. It vanishes when the displacement noraal to the boundary 
vanishes. 

We also found that dw FC cannot be ignored for displacements with 
t * n><0. Thus, there are two types of displacements? the isolated coronal 
■odes, with t • n * 0, for which 6w s and 6w pc can legitimately be set equal to 
zero, and the displacements with ? * n f 0. For the latter, the stability 
problem consists of jointly minimizing 5 W F q, 6Wp^ , and fiWg. This is not only 
impractical, in view of the lack of available equilibrium models, but is counter 
to the concept of an isolated coronal MHD instability. Virtually all papers on 
the subject have been concerned with isolated coronal instabilities, whether or 
not they treated the lower boundary condition in a way consistent with that 
idea. 

if the volume term alone is minimized, as was done by Hood (1984a,b) 

the result can be used to give a sufficient, but not necessary, condition for 
■tability of the isolated coronal modes. We showed that Hood's minimization 
of dw FC for a particular set of equilibria (Hood 1984b) led to nonvanishing 

t A -4- 

• n. Hie imposition of 5 • n ■= 0 on the boundary requires that Hood's trial 
functions A 1 have odd parity. We would argue, therefore, that some of the 
equilibria that Hood predicted are unstable are actually stable, according to 
the upper and lower^fluid model. 

The rigid boundary condition with t “ 0 has a vanishing surface term. 

•f 

Assuming the rigid boundary condition with £ * 0 leads to a self-consistent 
problem in which dw^ alone is minimized . This seems to be the simplest 
approach to treating the corona as an Isolated system. The full problem. 




t*v 


> V. 4. - 


• f. 


involving thermal exchange and dynamical forcing by motions of the fieldline 

* l 

endpoints, will have to be explored by other methods than the MHD energy 
principle. 
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